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ABSTRACT 

Experimental  testing  of  a  generic  submarine  model  has  been  conducted  in  the  DSTO  low  speed 
wind  tunnel.  These  tests  aimed  to  investigate  the  off -body  flow-field  around  a  generic  submarine 
shape  in  its  bare-hull  configuration  by  performing  smoke  flow  visualisation  and  gathering 
quantitative  velocity  data  using  particle  image  velocimtery  (PIV).  The  PIV  measurements  were 
compared  with  computational  fluid  dynamic  (CFD)  results  and  good  agreement  was  obtained.  The 
tests  also  provided  the  opportunity  to  refine  the  application  of  the  PIV  technique  in  a  large-scale 
wind  tunnel  facility.  These  experimental  data  will  complement  computational  and  experimental 
hydrodynamic  analysis  of  the  generic  submarine  shape. 
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Smoke  Flow  Visualisation  and 
Particle  Image  Velocimetry  Measurements  over  a 
Generic  Submarine  Model 


Executive  Summary 


The  Defence  Science  and  Technology  Organisation  (DSTO)  Undersea  Warfare  Corporate 
Enabling  Research  Program  (CERP)  is  aimed  at  improving  the  general  knowledge  and 
understanding  of  the  behaviour  of  complex  flows  around  modern  submarine  shapes.  The 
research  involves  the  use  of  high-fidelity  computational  fluid  dynamic  (CFD)  models,  as 
well  as  experimental  hydrodynamic  and  aerodynamic  test  techniques  to  investigate  a 
common  generic  shape.  The  extensive  dataset  gathered  from  the  variety  of  sources  will 
enhance  DSTO  knowledge  and  understanding  of  these  complex  flows  and  in  turn  serves 
as  validation  data  for  computational  methods. 

This  report  documents  the  wind  tunnel  testing  of  the  generic  submarine  model  in  the 
DSTO  low  speed  wind  tunnel  (LSWT)  in  early  2011.  The  aim  of  the  study  was  to 
characterise  the  flow-field  around  a  generic  submarine  shape  in  its  bare-hull  configuration, 
that  is,  without  the  hull-casing,  fin  and  control  surfaces  attached.  Experiments  were 
conducted  on  a  1.35  m  long  aluminium  model  with  a  diameter  of  0.185  m.  Smoke  flow 
visualisation  was  used  to  provide  qualitative  information  for  the  flow-field  as  the  model 
was  pitched  and  yawed.  Particle  Image  Velocimtery  (PIV)  then  provided  two-dimensional 
quantitative  velocity  measurements  along  the  centre-line  of  the  model  at  various 
incidences. 

Smoke  visualisation  results  revealed  that  the  flow-field  at  incidence  was  characterised  by 
the  formation  of  cross-flow  separation  near  the  tail,  which  migrated  forward  with 
increasing  incidence.  The  PIV  technique  provided  instantaneous  and  mean  flow-fields, 
from  which  maps  of  the  vorticity  and  Reynolds  stress  could  be  derived.  With  increasing 
incidence  the  PIV  derived  streamlines  along  the  upper  surface  of  the  model  displayed  a 
local  downwash  which  was  consistent  with  the  smoke  visualisation  results.  The  PIV 
results  were  compared  with  CFD  results.  The  streamwise  velocity  and  vorticity  showed 
good  agreement  with  the  CFD  results.  PIV  proved  to  be  a  useful  and  reliable  technique  for 
characterising  the  off-body  flow-field  around  the  generic  submarine  shape. 

Further  experimental  testing  of  the  model  in  the  DSTO  LSWT  is  planned,  including 
additional  PIV  measurements  with  the  hull-casing,  fin  and  the  control  surfaces  attached. 
The  implementation  of  a  range  of  improvements  suggested  in  this  report  will  provide 
DSTO  researchers  with  an  enhanced  PIV  capability  for  future  test  campaigns. 
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1.  Introduction 


Researchers  at  the  Defence  Science  and  Technology  Organisation  (DSTO)  have  used  smoke 
flow  visualisation  and  particle  image  velocimetry  (PIV)  to  assess  the  off-body  flow 
characteristics  around  a  model  of  a  generic  submarine.  This  work  was  conducted  in  the 
DSTO  low  speed  wind  tunnel  (LSWT),  with  the  secondary  aim  of  refining  the  application 
of  PIV  methods  when  applied  in  a  large-scale  wind  tunnel  facility. 

Smoke  flow  visualisation  was  used  to  provide  a  qualitative  picture  of  the  flow  around  the 
submarine  model,  and  provided  insight  into  the  nature  of  the  flow-field.  The  PIV 
measurements  were  then  targeted  at  areas  of  interest,  which  included  the  tail  and  wake 
region  of  the  model.  This  aft  region  is  typically  where  the  propeller  is  located,  and  at  full- 
scale  certain  inflow  conditions  can  lead  to  flow  induced  vibrations  and  noise. 

This  work  was  conducted  under  the  DSTO  Undersea  Warfare  -  Corporate  Enabling 
Research  Program  (CERP),  and  is  part  of  a  large  research  activity  to  improve  general 
knowledge  and  understanding  of  the  behaviour  of  complex  flows  around  modern 
submarine  shapes.  This  experiment  is  an  extension  of  the  force  and  moment  tests  of  the 
generic  submarine  model  conducted  in  the  DSTO  LSWT  [1].  In  addition  to  these 
experimental  tests,  the  same  model  will  also  be  tested  in  a  hydrodynamic  facility,  with  the 
results  compared  against  those  derived  using  computational  fluid  dynamics  (CFD) 
methods. 

This  report  documents  the  smoke  flow  visualisation  and  PIV  testing  in  the  LSWT 
including  descriptions  of  the  experimental  equipment  used,  the  wind  tunnel  model,  and 
the  test  facility.  The  test  methodology  and  the  test  conditions  are  also  presented,  together 
with  selected  results.  The  PIV  results  are  compared  qualitatively  with  the  smoke 
visualisation  and  quantitatively,  with  CFD  results.  A  number  of  key  recommendations  for 
future  testing  are  also  reported. 


2.  Experimental  Test  Facility  and  Model 

2.1  DSTO  Low  Speed  Wind  Tunnel 

Experiments  were  conducted  in  the  DSTO  LSWT  facility  located  at  Fishermans  Bend  in 
Melbourne.  The  facility  is  a  conventional,  continuous  flow,  closed-circuit,  single-return 
wind  tunnel  with  a  contraction  ratio  of  4:1.  The  test  section  has  an  irregular  octagonal 
shape  and  measures  2.74  m  (wide)  x  2.13  m  (high)  x  6.56  m  (long).  The  tunnel  is  powered 
by  a  660  kW  electric  motor  driving  a  3.96  m  diameter  8-bladed  fan.  The  facility  is  capable 
of  a  maximum  airspeed  of  approximately  100  m/s  when  the  test  section  is  empty. 
Honeycomb  consisting  of  triangle  cells  with  dimensions  48  mm  by  41  mm,  and  127  mm 
long,  is  installed  in  the  settling  chamber.  The  freestream  u-component  turbulence  intensity 
(u')  is  less  than  0.4%,  and  the  v  and  w  component  intensities  (v',  w')  are  less  than  0.7%  [2]. 
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Uncertainty  in  the  measurement  of  the  wind  speed  in  the  LSWT  test  section  is  discussed  in 
Appendix  A. 


2.2  Generic  Submarine  Model 

The  generic  submarine  model  was  designed  and  manufactured  for  testing  in  both  the 
DSTO  LSWT,  and,  with  minor  modification,  a  water  tunnel.  The  model  length  (L)  is  1.35 
m  and  model  diameter  (D)  is  0.185  m.  The  model  is  machined  from  aluminium  and 
comprises  a  cylindrical  centre-body  with  an  ellipsoid  nose  and  a  streamlined  tail  section. 
When  fully  assembled  the  model  has  fineness  ratio  (L/D)  of  approximately  7.3.  Nominally 
representative  of  a  modern  submarine  shape,  this  generic  test  article  has  no  full-scale 
equivalent. 

2.2.1  Model  Configuration 

In  this  study  the  bare-hull  configuration  was  tested  and  refers  to  the  hull-form  of  the 
generic  submarine  model  with  no  appendages  attached,  i.e.  the  casing,  fin  and  control 
surfaces  were  not  fitted,  as  shown  in  Figure  1. 

2.2.2  Support  System 

During  these  tests,  the  model  was  fitted  with  a  'dummy'  balance  and  connected  to  a  single 
vertical  support  pylon.  A  dummy  balance  was  used  in  order  to  protect  the  operational 
strain  gauge  balance  from  the  PIV  seed  particles.  Figure  1  shows  the  model  mounted  on  a 
pylon  within  the  working  section  of  the  LSWT.  The  model  pitch  attitude  was  controlled 
via  a  trailing  support-arm  and  vertical  control  link  mechanism.  The  model  pitch  was 
measured  by  an  inclinometer.  Jewel  Instruments  LCF-300  unit,  fitted  inside  the  model.  The 
vertical  support  pylon  was  shrouded  by  a  detachable  aerodynamic  fairing,  and  mounted 
on  a  rotating  turn-table  allowing  the  model  to  be  yawed  relative  to  the  freestream  flow. 
The  model  yaw  attitude  was  measured  by  the  turn-table  encoder.  The  inclinometer  pitch 
angle  bias  and  precision  uncertainty  is  ±0.029°  and  ±0.023°,  respectively  [1].  For  the  turn¬ 
table  encoder  the  yaw  angle  uncertainty  is  ±0.1°  [1]. 

2.2.3  Blockage  Ratio 

The  overall  frontal  area  blockage  ratio  for  the  model  at  zero-incidence  was  estimated  to  be 
2.1%  of  the  test  section  cross-sectional  area.  This  value  represents  the  sum  of  the  frontal 
area  of  the  model  (i.e.  0.5%)  and  the  pylon  fairing  (i.e.  1.6%). 
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Figure  1  -  Generic  submarine  model  mounted  in  the  LSWT 


2.2.4  Boundary  Layer  Transition 

For  all  the  experiments  a  boundary  layer  trip  was  employed.  A  3  mm  wide  strip  consisting 
of  carborundum  grit  was  attached  circumferentially  around  the  nose  of  the  model  at  5%  of 
the  body  length.  An  empirical  method,  described  in  reference  [3],  was  used  to  determine 
the  appropriate  grit  size  (i.e.  size  80,  or  an  average  particle  diameter  of  0.21  mm)  for  this 
test  programme.  During  the  force  and  moment  experiments  a  sound  test  (using  a 
stethoscope  traversed  manually  within  the  boundary  layer)  and  a  thermal  imaging 
technique  confirmed  that  this  arrangement  successfully  transitioned  the  boundary  layer  to 
a  turbulent  state  [1]. 


3.  Experimental  Method 

3.1  Smoke  Flow  Visualisation 

Smoke  flow  visualisation  provided  a  qualitative  assessment  of  the  flow  around  the 
submarine  model.  An  Aerotech  Smoke  Generator  supplying  Shell  Ondina  EL  oil  to  a 
heated  probe  was  used  to  generate  the  smoke.  The  probe  was  attached  to  the  end  of  a  long, 
slender  rod  that  was  used  to  position  the  smoke  stream  into  the  flow-field  around  the 
submarine  model.  The  smoke  flow  patterns  were  simultaneously  recorded  on  video  at 
1080i  resolution  using  two  Sony  HXR-NX5P  cameras  from  two  orthogonal  viewing 
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positions  directly  above  and  to  the  side  of  the  test  section.  A  Pelco  standard  definition 
observation  camera  located  inside  the  tunnel,  downstream  of  the  model,  also  provided 
video  recording.  The  observation  camera  and  two  Sony  cameras  were  synchronised  using 
a  camera  flash  and  matching  the  appropriate  frame  from  each  camera  during  post¬ 
processing. 

Two  Dedolight  Universal  Light  Heads  were  used  for  lighting.  The  smoke  stream  and 
model  were  illuminated  by  positioning  one  of  the  light  heads  inside  the  test  section 
downstream  of  the  model  with  the  light  directed  along  the  body  axis  of  the  model.  In 
order  to  minimise  shadows  from  the  model  and  surface  reflections  back  to  the  recording 
cameras,  a  second  light  head  was  positioned  outside  the  test  section,  upstream  of  the 
model  at  approximately  45°  to  the  oncoming  flow.  To  improve  the  contrast  of  the  smoke 
stream,  the  working  section  walls  were  covered  with  matte  black  paper. 


3.2  Particle  Image  Velocimetry 

PIV  is  a  non-intrusive  optical  velocity  measurement  technique  yielding  planar  whole  field 
quantitative  measurements  of  instantaneous  and  mean  velocity.  The  technique  is 
described  schematically  in  Figure  2.  The  PIV  technique  requires  the  flow-field  to  be  seeded 
with  small  particles  which  accurately  follow  the  flow.  A  dual-pulse  laser  is  used  to 
illuminate  a  thin  sheet  of  the  flow  with  a  short  time  delay.  At,  between  the  two  pulses.  The 
light  scattered  by  the  seed  particles  is  recorded  by  a  digital  camera,  in  this  case  a  camera 
with  a  frame-straddling  CCD  sensor,  positioned  normal  to  the  laser  sheet.  By  prescribing 
the  time  delay  between  pulses  (At)  and  measuring  the  particle  displacement  (Ax)  between 
the  two  recorded  images,  the  velocity  field  can  be  computed  as  Ax/ At.  In  practice  the 
particle  displacement  is  computed  by  dividing  the  recorded  image  into  sub-areas  called 
interrogation  windows  (IW).  The  in-plane  displacement  of  the  particles  in  each  IW  is 
obtained  using  statistical  image  cross-correlation.  A  velocity  vector  map  over  the  whole 
image  area  is  obtained  by  repeating  the  cross-correlation  and  velocity  calculation  for  each 
IW.  A  detailed  description  of  the  PIV  technique  is  given  in  [4]. 

In  this  study  the  2D  plane  of  interest  was  illuminated  using  a  New  Wave  Solo  PIV  200XT 
Nd:YAG  dual  laser-head.  This  laser  emits  two  independently  controlled  beams  at  a 
wavelength  of  532  nm,  a  maximum  pulse  energy  of  200  mj  and  a  5  ns  pulse  width.  The 
laser  was  mounted  above  the  wind  tunnel  test  section  on  a  two  degree  of  freedom 
translating  system  as  illustrated  in  Figure  3.  The  laser  traverse  was  designed  using 
MiniTec  components.  An  optical  train  comprising  a  cylindrical  and  spherical  lens  was 
used  to  generate  the  light  sheet.  The  laser  sheet  was  aligned  along  the  model  centre-line 
and  parallel  to  the  flow  direction,  as  shown  in  Figure  3.  In  order  to  minimise  laser 
reflection  towards  the  CCD  camera,  it  was  necessary  to  shift  the  laser  sheet  by  1  mm  in  the 
positive  z-direction  (away  from  the  camera).  On  the  surface  of  the  model  the  laser  sheet 
was  measured  to  be  1.5  mm  thick. 
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Figure  2  -  PIV  measurement  principle  (image  from  Dantec  Dynamics  [5]) 


The  PIV  image  acquisition  system  consisted  of  a  TSI  PowerView™  Plus  11  megapixel  CCD 
camera,  fitted  with  a  200  mm  Micro  Nikkor  f/4  lens.  The  CCD  camera  has  a  pixel 
resolution  of  4008  x  2672,  a  pixel  size  of  9  jam  x  9  jam  and  a  frame  rate  of  2.07  Hz  in  frame 
straddling  mode.  Data  acquisition  and  processing  of  the  instantaneous  flow-field  images 
was  achieved  using  the  TSI  INSIGHT™  3G  software.  Ensemble  averaging  of  the 
instantaneous  vectors  fields  and  further  post-processing  was  then  performed  using  a 
Matlab  program  developed  in-house.  The  camera  was  mounted  on  a  three  degree  of 
freedom  translating  system  and  was  positioned  orthogonal  to  the  laser  sheet,  as  shown  in 
Figure  3.  The  camera  traverse  comprised  of  Thorlabs  XT95  rail  series  and  custom  built 
components. 

The  PIV  seed  particles  comprised  of  olive  oil  droplets  generated  through  a  Laskin  nozzle 
(TSI  6-jet  oil  droplet  generator)  which  was  connected  to  a  compressed  air  supply  and 
remotely  controlled  using  a  solenoid  valve.  The  wind  tunnel  circuit  was  seeded  globally 
with  the  oil  aerosol  introduced  downstream  of  the  model.  The  particles  have  a  specific 
gravity  of  approximately  0.92,  a  mean  diameter  of  approximately  1  pm,  and  a  particle 
relaxation  time  which  was  calculated  to  be  approximately  2.4  ps. 
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3.2.1  Axes  System 

The  axes  system  used  for  the  PIV  measurements  is  defined  in  Figure  3,  with  the  origin 
located  at  the  tail  of  the  submarine  model.  The  local  coordinate  system  is  defined  with  the 
x-direction  in-line  with  the  positive  freestream  direction,  positive  y-direction  is  upwards 
and  the  z-direction  is  positive  starboard. 


Laser  Traverse 


Camera  Traverse 


Figure  3  -  PIV  arrangement  with  the  model  mounted  in  LSWT 


3.3  Derived  Quantities 

Large  quantities  of  velocity  field  data  generated  from  the  PIV  measurements  allow  a 
number  of  important  flow  quantities  to  be  calculated,  such  as  spatial  velocity  gradients, 
out-of-plane  vorticity  and  Reynolds  stresses.  The  velocity  data  is  two-dimensional  in 
nature  and  has  a  temporal  functional  dependence,  consequently  derived  quantities  also 
contain  a  spatial  dependence  in  two-dimensions  and  a  temporal  dependence.  A  detailed 
uncertainty  analysis  of  each  flow  quantity  is  described  in  Appendix  C. 
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3.3.1  Velocity 

In  the  current  investigation,  experiments  were  performed  in  the  x-y  plane  (see  Figure  3). 
From  the  in-plane  displacements,  AX  and  AY,  calculated  from  the  PIV  images,  the  fluid 
velocity,  u  and  v  at  a  point  in  the  flow,  x,  y  and  at  time,  t,  may  be  obtained  by. 


x  Y  A XxM 

u(x,y,t)  = - - - 

A  t 

and 

,  v  A  YxM 

v(  x,  y,  t)  = - > 


(3.1) 


where  At  is  the  time  separation  between  successive  images  pairs  and  M  is  the  image 
magnification.  Applying  Eq  3.1  for  each  IW  within  a  PIV  image  pair  results  in  an 
instantaneous  two-dimensional  velocity  field. 


3.3.2  Mean  Velocity 

Analysis  of  turbulent  flow  data,  or  any  random  flow  data,  requires  the  definition  of  an 
averaging  operator  for  the  experiment.  The  ideal  average  is  the  ensemble  average,  which  is 
defined  as  the  mean  over  an  infinite  collection  of  realisations  [6].  A  flow  field  is  said  to  be 
statistically  stationary  in  time  if  the  ensemble  averages  of  all  properties  are  independent  of 
time  [6].  For  the  current  PIV  experiments  the  flow  is  assumed  statistically  stationary  in 
time  and  the  camera  is  triggered  to  take  samples  of  the  velocity  field  periodically  in  time. 
The  time  between  successive  frame  pairs  (here  ~0.5  seconds)  is  also  assumed  large  enough 
so  that  successive  measurements  are  statistically  uncorrelated.  An  ensemble  average  is 
then  approximated  by  summing  over  a  large,  but  finite  collection  of  these  samples  (N).  It  is 
assumed  that  the  ensemble  average  in  this  experiment  is  equal  to  the  time  average  and 
consequently  the  ensemble  average  for  infinite  samples.  The  effect  of  the  number  of  PIV 
samples  used  to  estimate  the  ensemble  average  in  terms  of  the  sampling  uncertainty,  or 
uncertainty  in  the  mean,  is  discussed  in  Appendix  C. 


The  mean  velocity  at  a  point  in  the  flow,  x,  y,  calculated  from  N  velocity  realisations  is 
defined  as, 

1  n 

u(  X,  y)  =  (u(  X,  7,  t))  =  —  £  u(  X,  tn ),  (3.2) 

where  (.)  denotes  the  expected  value  (or  ensemble  average)  over  N  realisation.  The  mean 
velocity  components  calculated  from  the  PIV  measurements  are  therefore, 

u( x ,  y ),  v( x ,  y ).  (3.3) 


UNCLASSIFIED 


7 


UNCLASSIFIED 

DSTO-TR-2944 

3.3.3  Vorticity  and  Velocity  Gradients 

The  measurement  of  the  velocity  components  in  Eq.  (3.2)  allowed  the  calculation  of  the  out 
of  plane  vorticity  component. 


a>z(x,  y,t ) 


dv  dll 
dx  dy 


(3.4) 


The  vorticity  was  derived  by  calculating  the  spatial  velocity  derivatives  (£^  ^H)  using  a 


dx  dy 


second-order  central  finite  difference  scheme  at  the  interior  points, 

df  _/(/-!)-/(/  +  !) 


dx  2A 

where /is  generic  function  and  i  is  the  spatial  coordinate  index  in  the  x  direction. 


(3.5) 


The  second  order  forward  difference. 


8/  -3/(Q  +  4/Q-  +  !)-/(»  +2) 

dx  2A 

and  backward  difference, 

df  f(i-2)-Af(i-l)  +  3f(i)  _ 
dx  2A 


(3.6) 


(3.7) 


schemes  were  used  at  the  boundary  of  the  region  of  interest.  For  interior  points  coinciding 
with  spurious  data,  second-order  finite  difference  could  not  be  used,  and  instead  a  second- 
order  forward  or  backward  difference  was  used.  Spurious  data  is  discussed  in  Section  3.4. 
Etabari  and  Vlachos  [7]  indicate  that  the  second-order  central  finite  difference  scheme 
gives  the  best  trade-off  in  performance  between  the  associated  bias  and  precision 
uncertainty  resulting  from  a  particular  computational  method  used  to  calculate  the 
velocity  derivatives. 

3.3.4  Mean  Vorticity 

The  mean  vorticity  at  a  point  in  the  flow,  x,  y,  calculated  from  N  velocity  realisations  is 
defined  as, 

<»z(x,y)  =  ((oz(x,y,t)).  (3.8) 

3.3.5  Turbulent  Quantities 

The  velocity  fluctuating  components  were  determined  as  follows, 

u'  (  x,y,t )  =  u(  x,y,t )  -u(  x,y  ) , 

v'(  x,y,t )  =  v(  x,y,t )  —  v(  x,y  )  • 


(3.9) 

(3.10) 
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The  in-plane  Reynolds  stress  at  a  point  in  the  flow,  x,  y,  calculated  from  N  velocity 
realisations  is  defined  as. 


(x,y)  =  (iu( x,y,t)-u(x,y )ju( x,y,t)-u(x,y )\  • 


(3.11) 


The  other  Reynolds  stress  components  (  v'2(x,y)  and  u'v'  (. x,y ))  can  similarly  be  obtained. 
The  turbulence  strength  or  root-mean-square  of  the  velocity  fluctuation  is 


uRMs(x’y)  =  ^u’  (x>y) 


(3.12) 


and  the  turbulence  intensity. 


URMs( x>  y )  /  U„  • 


(3.13) 


Similarly  the  v-component  turbulence  strength  and  turbulence  intensity  are  defined  as, 

(3.14) 


vrms( x>y)=\  v'2 (x,y),  and  vRMS( x, y)/Uao . 


3.4  Data  Processing  and  Validation  Criteria 

The  particle  displacements  are  derived  by  processing  the  PIV  image  pairs  using  the  Insight 
3G  software.  Within  each  IW  the  particle  displacement  are  computed  using  a  two-pass 
recursive  algorithm  with  an  IW  size  for  both  the  first  and  second  pass  of  64  x  64  pixels  and 
a  50%  grid  spacing  overlap.  Each  vector  in  the  resultant  vector  field  was  then  validated 
and  filtered  to  eliminate  spurious  (or  invalid)  vectors  which  may  be  the  result  of  low 
particle  seeding,  strong  velocity  gradients,  particles  moving  out-of-plane  due  to  three- 
dimensional  motion  or  a  noisy  image.  This  comprised  a  global  velocity  range  filter,  a  mean 
filter  using  a  neighbourhood  size  of  3  x  3  pixels  and  a  displacement  tolerance  of  2  pixels. 
Vectors  identified  as  invalid  were  then  replaced  with  interpolated  points  using  bi-linear 
interpolation. 

Liang  et  al.  [8]  state  that  for  well-optimised  PIV  experiments,  the  spurious  vector  count 
should  be  less  than  5%  of  the  total  number  of  measured  vectors.  In  this  study  if  a  single 
velocity  vector  field  contained  spurious  vector  numbers  greater  than  5%  of  the  total 
number  of  measured  vectors  it  was  not  included  in  the  ensemble  average.  This  resulted  in 
some  of  the  mean  flow-fields  being  averaged  over  less  than  the  number  of  acquired 
images.  The  final  number  of  samples  (N)  used  in  the  ensemble  average  for  each  test  case  is 
listed  in  Table  2.  Whilst  some  test  cases  had  fewer  samples  than  1000,  in  Section  4.4.2,  a 
sample  size  of  500  is  shown  to  be  sufficient  for  the  convergence  of  mean  flow  and 
turbulent  parameters. 
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In  some  cases  the  PIV  measurement  fields  contained  regions  of  invalid  data  located  in  the 
same  region  for  each  field.  Where  possible  these  regions  were  masked  out  leaving  only  the 
valid  PIV  data.  Further  discussion  regarding  the  factors  contributing  to  invalid  data 
regions  are  presented  in  Section  4.5.  After  post  processing  the  PIV  data,  three  test  cases 
were  discarded  as  they  were  found  to  have  a  significant  amount  (>10%)  of  spurious 
vectors  present  in  each  field.  Of  the  remaining  valid  fields,  statistical  analysis  revealed  that 
the  mean  number  of  spurious  vectors  in  a  single  processed  PIV  image  was  less  than  1.6% 
and  95%  of  the  fields  contained  less  than  2.8%  spurious  vectors.  Thus  for  the  majority  of 
test  cases,  the  PIV  measurements  were  of  a  high  quality. 


3.5  Test  Conditions 

The  generic  submarine  model  was  tested  at  various  angles-of-attack  (a)  and  angles-of- 
sideslip  (|3).  The  freestream  airspeed  (Uoo)  in  the  wind  tunnel  was  dictated  by  the 
experimental  technique  used.  The  smoke  flow  visualisation  technique  required  lower 
speeds  so  that  the  smoke  stream  did  not  disperse  too  quickly,  whereas  for  the  PIV 
measurement  higher  speeds  could  be  attained. 

The  smoke  flow  visualisation  tests  were  conducted  at  Uoo  =  5  m/s  and  the  following  model 
orientations: 

•  Pitch,  a  =  0°,  ±2.5°,  ±5°,  ±10°  and  ±15°  (at  zero  yaw) 

•  Yaw,  (3  =  0°,  5°,  10°,  20°  and  30°  (at  zero  pitch) 

•  Combined  pitch  and  yaw  of  -  15°  pitch/15°  yaw  and  15°  pitch/30o  yaw. 

The  test  matrix  for  the  PIV  tests  is  shown  in  Table  1.  The  majority  of  the  PIV 
measurements  were  conducted  at  Uoo  =  30  m/ s.  The  effect  of  Reynolds  number  based  on 
the  model  length  (ReL  =pUo0L/]i ),  was  studied  by  performing  PIV  measurements  at  Uoo  = 

20  m/ s,  30  m/s  and  40  m/ s  for  the  model  at  a  =  0°  and  |3  =  0°. 

Table  1  -  PIV  test  matrix 


an 

1100(111/  s) 

ReL 

0 

20 

1.8xl06 

2.5,  0,  -2.5,  5, 10, 15 

30 

2.7xl06 

0 

40 

3.6xl06 

The  PIV  test  conditions  and  associated  PIV  parameter  settings  are  defined  in  Table  2.  The 
field-of-view  column  refers  to  the  discrete  locations  on  the  model  where  measurements 
were  taken  and  are  detailed  in  Section  3.6.  The  time  intervals  between  laser  pulses  (At) 
were  based  on  a  requirement  that  the  maximum  particle  displacements  be  less  than  25%  of 
the  64  x  64  pixel  IW.  Three  test  cases  were  discarded  due  to  large  percentage  of  spurious 
vectors  (see  Section  3.4).  A  summary  of  the  PIV  measurement  settings  is  provided  in 
Appendix  B. 
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Table  2  -  PIV  test  conditions  and  flow  measurement  parameters  at  different  pitch  angles  at  zero 
yaw 


Field-of-view 

a(°) 

U co  (m/  s) 

At  ( ]is ) 

Number  of  fields 
(N) 

i 

-2.5 

29.9 

10 

1000 

2 

-2.5 

29.9 

10 

1000 

3 

-2.5 

29.9 

10 

1000 

1 

0 

19.8 

15 

960 

1 

0 

40.5 

6 

906 

1 

0 

30.4 

10 

827 

2 

0 

30.2 

10 

1000 

3 

0 

30.0 

10 

discarded 

1 

2.5 

30.3 

10 

910 

2 

2.5 

30.3 

10 

1000 

3 

2.5 

29.9 

10 

1000 

1 

5 

29.9 

10 

1000 

2 

5 

29.9 

10 

discarded 

3 

5 

29.9 

10 

1000 

1 

10 

29.9 

10 

1000 

2 

10 

29.9 

10 

999 

3 

10 

29.9 

10 

1000 

1 

15 

30.0 

10 

986 

2 

15 

30.3 

13 

discarded 

3 

15 

30.3 

13 

1000 

3.6  Particle  Image  Velocimetry  Measurement  Locations 

The  image  field-of-view  is  the  visible  extent  of  the  flow  field  captured  on  the  camera 
sensor.  It  is  desirable  for  the  field-of-view  to  be  sufficiently  large  so  as  to  capture  the  entire 
region  of  interest  in  a  single  realisation.  Conversely,  high  spatial  resolution  is  required  in 
order  to  adequately  resolve  flow  features  around  the  model.  These  requirements  are 
conflicting,  since  a  high  spatial  resolution  results  in  a  smaller  field-of-view. 

In  addition  to  the  above  constraints,  another  important  PIV  requirement  is  to  avoid  peak¬ 
locking  (also  known  as  pixel-locking).  Peak-locking  is  the  tendency  of  PIV  measurements 
to  bias  towards  integer  displacements.  This  occurs  when  the  particle  images  are  under¬ 
resolved  and  the  peak  fit  applied  to  the  cross-correlation  function  (a  three-point  Gaussian) 
is  unable  to  accurately  determine  the  correlation  peak  to  sub-pixel  accuracy  [9].  For  a 
three-point  Gaussian  peak  fit  applied  to  the  cross-correlation  function,  the  optimum 
particle  image  diameter  for  PIV  measurements  is  between  2  to  3  pixels  [4].  Practically  in 
this  test,  such  a  particle  image  diameter  was  not  possible  due  to  the  following  factors: 

•  the  diameter  of  the  olive  oil  droplets  were  small  (nominally  1  jam), 

•  the  image  magnification  was  limited  as  the  camera  was  located  relatively  far  from 
the  object  measurement  plane,  and 

•  the  maximum  focal  length  lens  available  was  a  200  mm  Micro-Nikkor. 
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A  typical  approach  to  improve  the  particle  image  size  is  to  increase  the  camera  lens  f-stop, 
which  increases  the  diffraction  limited  particle  diameter.  The  drawback  with  this  approach 
is  that,  in  order  to  maintain  the  same  exposure,  each  increase  in  f-stop  requires  a  doubling 
of  laser  light  intensity.  Alternatively,  to  achieve  a  larger  image  particle  diameter  (and 
higher  spatial  resolution)  an  imaging  lens  with  a  longer  focal  length  could  be  used,  but 
was  not  available  for  this  experiment.  Peak-locking  for  the  current  experiment  is  assessed 
in  Section  4.4.3. 

For  this  study,  in  order  to  image  an  overall  larger  area  and  have  sufficient  resolution  to 
avoid  peak-locking,  three  high  resolution  field-of-views  were  blended  together.  The 
distance  from  the  centre  of  the  test  section  to  the  PIV  camera  was  approximately  1.6  m  and 
using  the  200  mm  lens,  a  reproduction  ratio  (i.e.  ratio  of  subject  size  on  sensor  plane  to  the 
actual  subject  size)  of  6.67  resulted.  The  corresponding  image  spatial  resolution  was 
60.03  pm/ pixel  and  the  field-of-view  for  each  image  was  240  mm  x  160  mm. 

The  PIV  measurements  were  focused  around  the  tail  and  wake  region  of  the  generic 
submarine  model.  The  locations  of  the  three  field-of-views  over  which  measurements  were 
planned  are  shown  in  Figure  4.  However,  due  to  the  optical  distortions  produced  by  the 
test  section  window  (described  in  Section  4.5.2),  these  locations  could  not  always  be 
obtained  for  each  measurement  case.  Overlapping  the  field-of-views  allowed  the  data  to 
be  blended  together  during  post-processing  using  the  approach  described  in  Kumar  et.  al 
[10].  The  amalgamation  of  the  three  field-of-views  is  shown  in  Figure  5  with  respect  to  the 
entire  model  body.  Blending  multiple  high  resolution  fields  is  common  practise  in  studies 
using  PIV,  and  has  been  reported  in  submarine  research  [11]. 
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4.  Results 


4.1  Smoke  Flow  Visualisation 

Smoke  visualisation  provided  qualitative  insight  into  the  flow  around  the  submarine 
model  and  assisted  in  the  interpretation  of  the  PIV  measurements.  Figure  6  shows  smoke 
flow  visualisation  patterns,  on  the  starboard  side  of  the  model,  with  increasing  a.  The  inset 
image  shows  the  camera  view  from  above  the  test  section  taken  at  the  same  instance  in 
time.  At  a  =  0°,  Figure  6(a)  shows  the  flow  is  everywhere  attached  to  the  model.  At  a  = 
2.5°,  the  flow  around  the  side  of  the  hull  remains  attached  and  exhibits  a  small  upward 
curvature.  As  a  is  increased  further,  the  curvature  of  the  flow  increases  and  begins  further 
upstream.  At  a  =  5°,  there  is  evidence  of  separation  at  the  tail.  With  increasing  a,  this 
cross-stream  separation  migrates  upstream.  The  resultant  separated  flow  rolls  up  into  a 
longitudinal  vortex  on  each  side  of  the  body  (as  depicted  schematically  for  a  prolate 
spheroid  in  Figure  7).  These  longitudinal  vortices  are  stronger  and  more  pronounced  as 
more  smoke  is  entrained  towards  the  core  centre  and  the  location  where  they  detach  from 
the  body  moves  further  upstream  with  increasing  a. 

The  flow  patterns  observed  in  the  LSWT  flow  visualisation  study  of  the  hull  body  are 
consistent  with  studies  of  prolate  spheroids  of  various  axis  ratios  [12-16].  The  flow  around 
a  prolate  spheroid  is  characterised  by  cross-stream  separation  that  forms  on  the  tail  at  low 
a  and  then  migrates  forward  and  leeward  with  increasing  a.  The  inward  roll  up  of  the 
longitudinal  vortex,  induces  a  downward  velocity  component  along  the  vertical  symmetry 
plane  of  the  model  (refer  to  Figure  7).  As  the  model  was  pitched,  this  downward  velocity 
component  resulted  in  a  centre-line  flow  which  remained  attached  for  all  a  tested.  An 
example  of  the  attached  flow  along  the  vertical  symmetry  plane  of  the  model  is  shown  in 
Figure  8,  for  a  =  15°.  The  two  images  are  the  same  instance  in  time  taken  from  the 
starboard  side  video  camera  and  the  observation  camera  inside  the  LSWT. 

The  generic  submarine  model  in  the  bare  hull  configuration  is  axisymmetric,  and  in  the 
absence  of  support  interference  effects,  (3  is  analogous  to  a.  This  was  evident  in  the  smoke 
flow  visualisation  with  the  large-scale  flow  behaviour  as  p  is  increased.  However,  unlike 
pitch  motion,  increasing  p  showed  the  longitudinal  vortices  were  present  at  an  oblique 
angle.  Due  to  the  mounting  arrangement,  higher  angles  of  p  (up  to  30°)  compared  to  a 
could  be  achieved.  At  p  =  30°  for  example.  Figure  9  show  the  smoke  to  be  entrained  in  two 
distinct  vortices  which  are  approximately  symmetric  about  the  horizontal  plane. 

The  smoke  flow  visualisation  results  clearly  show  the  influence  of  the  mounting 
arrangement,  particularly  on  the  underside  of  the  model,  toward  the  aft  end.  A  turbulent 
wake  downstream  of  the  pylon-body  junction  was  apparent,  as  well  as  flow  separation 
and  in  some  instances,  vortex  shedding.  With  the  model  at  a  =  0°  and  the  smoke  probe 
positioned  in  the  aft  fairing  cut-out,  a  recirculation  region  was  also  evident  and  is  shown 
in  Figure  10.  The  flow  around  the  fairing  separates  at  the  cut-out,  and  a  portion  of  the  flow 
enters  the  recirculation  region.  On  exiting  the  recirculation  region,  the  flow  is  disturbed  by 
the  pitch  actuator  arm.  The  flow  interference  effect  appears  to  be  greatest  as  the  model  is 
yawed.  For  large  side-slip  angles  the  pitch  actuator  arm  is  no  longer  in  the  pylon  and 
fairing  wake  and  becomes  exposed  to  the  oncoming  flow.  An  example  of  vortex  shedding 
from  the  pitch  actuator  arm,  at  p  =  30°,  is  shown  in  Figure  11.  An  unsteady  flow-field 
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downstream  of  the  vertical  support  pylon-model  body  junction  was  also  observed  in  the 
wool  tuft  flow  visualisation  experiments  performed  by  Quick  et  al.  [1]. 

(a) 


(b) 


(c) 


(d) 


Figure  6  -  Smoke  flow  visualisation  patterns  on  the  starboard  side  of  the  submarine  model  with 
increasing  angle-of-attack  (a)  a  =  0°,  (b)  a  =  2.5°  (c)  a  =  5°  and  (d)  a  =  15° 

1Inset  image  show  the  view  from  above  the  model  at  the  same  instant  in  time. 
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Figure  7  -  Cross-flow  separation  on  a  prolate  spheroid  at  pitch  (from  Chesnakas  and  Simpson  [16]) 


Figure  8  -  Smoke  flow  visualisation  at  a  =  15°  and  f  =  0°,  showing  attached  flow  along  the  centre¬ 
line  of  the  model 


Figure  9  -  Smoke  flow  visualisation  at  a  =  0°  and  f  =  30°,  showing  roll-up  of  the  vortex  sheets  into 
two  distinct  vortices 
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Figure  10  -  Smoke  flow  visualisation  at  a  =  0°  and  (3  =  0°,  with  probe  positioned  at  the  aft  fairing 
cut-out  showing  a  recirculation  region  and  an  unsteady  wake 


Figure  11  -  Smoke  flow  visualisation  at  a  =  0°  and  (3  =  30°,  showing  vortex  shedding  from  the 
pitch  actuator  arm 


4.2  PIV  Mean  Flow-field 

The  mean  PIV  flow-field  results  were  obtained  by  an  ensemble  average  of  the 
instantaneous  vector  fields  over  the  number  of  samples  taken  for  each  case,  N,  as 
summarised  in  Table  2  in  Section  3.5. 

4.2.1  Flow-field  at  Low  Angles-of- Attack 

For  a  =  0°,  Figure  12(a)  shows  the  non-dimensional  streamwise  velocity  field  (u/U^), 
together  with  the  streamline  pattern.  Attached  flow  passes  over  the  tail  surface  of  the 
submarine  model  with  the  streamlines  showing  a  small  downward  velocity  component. 
The  out-of-plane  non-dimensional  vorticity  (& ZL/Uoo)  is  presented  in  Figure  12(b).  The 
boundary  layer  vorticity  expands  perpendicular  to  the  model  surface  with  downstream 
distance  and  forms  a  free  shear  layer  beyond  the  tail.  The  region  of  higher  vorticity 
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extends  further  downstream  with  decreasing  a.  The  Reynolds  stress  (-u'v'/U^)  contours 
are  presented  in  Figure  12(c).  The  Reynolds  stress  plot  displays  the  level  of  turbulence 
present  in  the  flow,  and  shows  strong  intensities  at  the  upstream  end  of  the  field-of-view 
and  in  the  free  wake  region.  It  should  be  noted  that  in  this  and  all  other  Reynolds  stress 
contour  plots  presented,  the  magnitude  of  the  Reynolds  stress  fluctuations  are  low 
compared  to  the  experimental  uncertainty  (see  Appendix  C).  The  Reynolds  stress  contour 
plots  therefore  only  provide  a  qualitative  assessment  of  the  flow.  Strategies  on  how  to 
reduce  the  PIV  measurement  uncertainty  are  described  in  Section  4.4.1. 

The  PIV  results  at  a  =  ±2.5°  display  similar  features  to  a  =  0°  and  are  shown  in  Figure  13 
and  Figure  14.  In  both  cases,  the  boundary  layer  over  the  tail  surface  remains  attached, 
with  the  stagnation  point  located  at  the  model  trailing  edge.  At  -2.5°  the  vorticity  over  the 
top  of  the  model  displays  greater  intensity,  while  at  +2.5°  it  is  more  diffuse.  In  the  absence 
of  model  support  interference  effects,  the  flow  features  at  a  of  +2.5°  and  -2.5°  should  be  a 
mirrored  image  about  the  y/L  =  0  line.  However,  both  exhibit  increased  levels  of 
turbulence  in  the  lower  wake  region  due  to  the  flow  unsteadiness  generated  by  the  vertical 
support  pylon  as  identified  in  the  flow  visualisation. 

4.2.2  Flow-field  at  High  Angles-of- Attack 

The  PIV  results  as  the  model  is  pitched  at  5°,  10°  and  15°  are  shown  in  Figure  15,  Figure  16 
and  Figure  17,  respectively.  With  increasing  a,  the  streamlines  over  the  top  of  the  model 
are  vectored  down  along  the  surface.  This  is  consistent  with  the  flow  visualisation  results 
at  these  a,  two  longitudinal  vortices  roll  up  inward  towards  the  surface  of  the  model  and 
induce  a  downward  velocity  along  the  vertical  symmetry  plane  of  the  model.  This 
downwash  is  also  depicted  in  Figure  7  and  is  shown  by  the  half-saddle  in  cross-section. 
For  the  tested  a,  this  downward  velocity  ensures  the  flow  remains  attached  along  the 
vertical  symmetry  plane  of  the  model. 

The  vorticity  contours  show  that,  as  a  increases  from  5°  to  10°,  the  boundary  layer 
vorticity  on  the  upper  surface  extends  further  downstream.  With  increasing  a  the  shear 
layer  vorticity  is  " compressed"  closer  to  the  boundary  of  the  model.  At  a  =  10°  and  15°,  a 
thin  streak  of  positive  vorticity  is  also  evident  above  the  model,  displaced  from  the  model 
surface  at  ~y/L  =  0.08.  Since  the  flow  visualisation  at  these  a  confirmed  the  presence  of 
cross-stream  separation  and  a  flow  topology  similar  to  the  sketch  in  Figure  7,  this  streak 
may  be  an  artefact  produced  by  the  local  junction  flow  between  the  counter-rotating 
vortices  (i.e.  the  off -body  node  line)  near  the  model  symmetry  plane. 

The  Reynolds  stress  contour  plots  show  streamwise  structures  which  become  more 
concentrated  with  increasing  a. 
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Figure  12  -  (a)  Streamline  topology  overlayed  with  non-dimensional  streamwise  velocity ,  (b)  non- 
dimensional  out-of-plane  vorticity,  (c)  Reynolds  stress  contours  at  a  =  0° 
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Figure  13  -  (a)  Streamline  topology  overlayed  with  non-dimensional  streamwise  velocity ,  (b)  non- 
dimensional  out-of-plane  vorticity,  (c)  Reynolds  stress  contours  at  a  =  -2.5° 
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Figure  14  -  (a)  Streamline  topology  overlayed  with  non-dimensional  streamwise  velocity,  (b)  non- 
dimensional  out-of-plane  vorticity ,  (c)  Reynolds  stress  contours  at  a  =  2.5° 


UNCLASSIFIED 


21 


DSTO-TR-2944 


UNCLASSIFIED 


Figure  15  -  (a)  Streamline  topology  overlayed  with  non-dimensional  streamwise  velocity,  (b)  non- 
dimensional  out-of-plane  vorticity,  (c)  Reynolds  stress  contours  at  a  =  5° 
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Figure  16  -  (a)  Streamline  topology  overlayed  with  non-dimensional  streamwise  velocity,  (b)  non- 
dimensional  out-of-plane  vorticity,  (c)  Reynolds  stress  contours  at  a  =  10° 
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Figure  17  -  (a)  Streamline  topology  overlay ed  with  non-dimensional  streamwise  velocity,  (b)  non- 
dimensional  out-of-plane  vorticity,  (c)  Reynolds  stress  contours  at  a  =  15° 
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4.2.3  Reynolds  Number  Sensitivity  Study 

A  Reynolds  number  sensitivity  study  was  conducted  at  a  =  0°  and  for  field-of-view  1  (see 
Figure  4).  Freestream  velocities  of  20  m/s,  30  m/s  and  40  m/s  were  tested,  which 
correspond  to  Reynolds  numbers,  based  on  the  model  length,  of  1.8xl06,  2.7xl06  and 
3.6xl06,  respectively.  The  non-dimensional  streamwise  velocity  contours,  the  out-of-plane 
non-dimensional  vorticity  and  the  Reynolds  stress  contours  are  presented  for  each 
Reynolds  number  in  Figure  18.  The  flow  field  appears  independent  of  ReL  for  the  range  of 
velocities  tested,  as  the  contour  plots  for  all  three  cases  display  similar  flow  patterns.  An 
exception  is  the  Reynolds  stress  plot  for  Uoo  =  40  m/s,  where  two  bands  of  concentrated 
negative  Reynolds  stress  appear  near  the  surface  of  the  model  and  then  dissipate  with 
downstream  distance.  These  low  concentration  bands  are  also  evident  in  the  vorticity 
contour  plots  and  are  more  intense  with  increasing  ReL.  It  remains  unclear  if  these  bands 
are  physical,  due  to  the  relatively  high  measurement  uncertainty  for  these  quantities.  As 
discussed  in  Section  4.4.1,  in  the  uncertainty  in  PIV  measurements  can  be  reduced  by 
increasing  the  spatial  resolution  and  increasing  the  sample  size  used  in  the  ensemble 
average. 

For  this  study  there  is  a  disparity  between  the  Reynolds  numbers  in  the  smoke  flow 
visualisation  and  PIV  tests.  The  smoke  visualisation  testing  was  conducted  at  ReL  =  0.45  x 
106,  whereas  the  PIV  experiments  were  predominantly  conducted  at  ReL  =  2.7  x  106.  In  the 
case  of  flow  over  a  prolate  spheroid  of  axis  ratio  6:1,  it  was  discovered  by  Ahn  [17]  and 
then  confirmed  by  Wetzel  and  Simpson  [18],  that  a  critical  Reynolds  number  exists  above 
which  the  flow-field  is  independent  of  the  Reynolds  number.  They  found  that  the  cross- 
flow  separation  was  insensitive  to  Reynolds  number  above  ReL  =  2.5  x  106.  For  the  current 
study  it  is  unclear  if  such  a  critical  Reynolds  number  exists.  The  Reynolds  number  effect 
study  was  conducted  at  a  =  0°  and  not  at  an  incidence  where  the  cross-stream  separation 
is  likely  to  be  evident  (i.e.  a  >  5°).  Further  investigation  is  required  in  order  to  assess  if  a 
critical  Reynolds  number  exists  for  the  generic  submarine  shape  tested  here.  Such 
investigations  should  utilise  surface  flow  visualisation  techniques  as  well  as  other  test 
techniques  which  will  enable  the  measurement  of  the  cross-stream  flows. 
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Figure  18  -  Ensemble-averaged,  (a)  non-dimensional  streamwise  velocity  (u/Lly),  (b)  non- 
dimensional  out-of-plane  vorticity  (a)zL/Ux) ,  (c)  Reynolds  stress  (-u'v'/Ul) 
contours  for  field-of-view  1,  a  =  0°  and  for  freestream  velocities  of 20,  30  and  40  m/s. 


4.3  Comparison  of  Experimental  and  Computational  Results 

As  part  of  the  verification  and  validation  of  the  PIV  data,  the  experimental  results  were 
compared  with  available  computational  fluid  dynamics  (CFD)  results  [19].  The  CFD  model 
of  the  generic  submarine  including  the  support  fairing  and  pitch  actuator  arm  is  shown  in 
Figure  19.  To  reduce  the  complexity  of  the  CFD  model  the  wind  tunnel  test  section,  the 
central  vertical  support  pylon  connecting  the  support  fairing  with  the  model,  and  the  pitch 
strut  vertical  linkage  were  not  modelled.  Whilst  the  CFD  simulation  is  a  simplification  of 
the  experimental  setup,  it  was  shown  in  Quick  et  al.  [1]  that  the  forces  predicted  correlate 
well  with  the  experimental  results. 

In  the  CFD  model  an  unstructured  tetrahedral  mesh  was  used,  with  the  addition  of  a 
prism  layer  near  the  surface  to  capture  the  boundary  layer.  An  initial  cell  height  of  0.001  m 
coupled  with  a  radial  growth  rate  of  1.15  using  20  prism  levels  was  deemed  sufficient  to 
capture  the  boundary  layer  development  [20].  The  average  surface  cell  i/+  value  was 

approximately  80  and  the  mesh  size  was  approximately  lxlO6  cells.  A  standard  k-co 
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turbulence  model  was  used  with  a  wall  function  approach.  The  computational  analysis 
assumed  a  fully  turbulent  boundary  layer. 

For  the  model  at  a  =  10°  a  comparison  between  CFD  and  PIV  derived  non-dimensional 
streamwise  velocity,  vorticity  and  Reynolds  stress  are  made  in  Figure  20,  Figure  21  and 
Figure  22,  respectively.  For  the  non-dimensional  streamwise  velocity  a  good  qualitative 
agreement  is  evident  with  both  results  displaying  similar  flow  features.  However,  the 
wake  in  the  PIV  plots  is  broader  and  extends  further  downstream.  Differences  in  the 
vorticity  and  Reynolds  stress  are  more  notable,  with  the  PIV  measurements  indicating 
larger  concentrated  regions  and  higher  turbulence  in  the  wake  of  the  model.  Discrepancies 
may  be  attributed  to  inaccuracies  and  uncertainties  in  both  methods.  One  notable 
discrepancy  is  the  central  vertical  support  pylon,  which  was  not  modelled  in  the  CFD  and 
may  account  for  some  of  the  differences  observed  between  the  two  results. 


Figure  19  -  The  CFD  surface  model  of  the  generic  submarine  including  the  support  fairing  and 
actuator  arm  [19] 

To  provide  a  quantitative  comparison  between  the  PIV  and  CFD  results,  profiles  were 
extracted  at  x/L  =  -0.059,  -0.02,  0.011  and  0.05  relative  to  the  tail  tip  for  the  model  at  a  =  0°. 
Figure  23  presents  plots  of  the  mean  non-dimensional  streamwise  velocity  at  these  four 
locations  for  the  model  at  a  =  0°.  The  error  bands  represent  the  uncertainty  in  the  mean 
flow  PIV  measurements  (refer  to  Appendix  C  for  the  uncertainty  analysis).  Overall  the 
streamwise  velocity  profiles  compare  well.  The  streamwise  velocity  predict  by  CFD  is  less 
than  that  measured  using  PIV,  except  for  x/L  =  0.05  which  shows  an  underestimate  of  the 
turbulence  quantities  in  the  wake  region.  For  x/L  =  -0.059  and  -0.02  the  largest  difference 
is  observed  near  the  wall,  where  the  velocity  gradients  and  therefore  the  PIV  bias 
uncertainties  are  greatest. 

Figure  24  presents  the  mean  non-dimensional  out-of-plane  vorticity  at  the  same  four 
locations.  The  oscillatory  nature  of  the  PIV  profiles  is  not  evident  in  the  CFD  results, 
however  the  CFD  profiles  do  fall  within  the  experimental  error  bands  of  the  PIV  results. 
The  oscillatory  behaviour  is  likely  due  to  the  relatively  large  uncertainty  in  the  PIV 
vorticity  measurements  (see  Appendix  C).  Further  investigation  is  required  to  ascertain  if 
this  behaviour  is  physical. 
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Figure  21  -  Comparison  between  CFD  generated  non-dimensional  out-of-plane  vorticity  (FCzl/u„) 
and  PIV  results ,  at  a  =  10° 
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Figure  22  -  Comparison  between  CFD  generated  Reynolds  stress  (-u'v'/Ul)  contours  and  PIV 
results,  at  a  =  10° 
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Figure  23  -  CFD  and  PIV  non-dimensional  streamwise  velocity  profiles  at  a  =  0°,  at  locations,  (a) 
x/L  =  -0.059,  (b)  x/L  =  -0.02,  (c)  x/L  =  0.011  and  (d)  x/L  =  0.05 
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Figure  24  -  CFD  and  PIV  non-dimensional  out-of-plane  vorticity  profiles  at  a  =  0°,  at  locations  (a) 
x/L  =  -0.059,  (b)  x/L  =  -0.02,  (c)  x/L  =  0.011  and  (d)  x/L  =  0.05 
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4.4  Assessment  of  Data  Quality 

4.4.1  Measurement  Uncertainties 

The  measurement  uncertainties  in  the  instantaneous,  mean  and  fluctuating  velocity 
components  were  calculated  using  the  approach  described  in  [21]  and  are  listed  in  Table  3. 
The  uncertainties  in  the  mean  and  fluctuating  velocity  components  are  based  on  random 
variable  analysis  for  a  95%  confidence  interval.  The  vertical  velocity  component  (v)  is 
assumed  to  have  identical  uncertainty  as  the  streamwise  velocity  component  (u).  The 
uncertainties  presented  are  of  the  same  order  as  those  in  other  published  PIV  data  [22].  A 
detailed  PIV  uncertainty  analysis  is  presented  in  Appendix  C. 

Table  3  -  Uncertainty  of  the  normalised  instantaneous  and  ensemble-averaged  velocity ,  out-of- 
plane  vorticity  and  Reynolds  stress 


Ucc 

(m/s) 

N 

Uncertainty 

Eu 

£4 

ejl 

Ux 

E- 

co  ZL 

Ux 

E — 

u  V 

uf 

30 

1000 

±7.4xl0'3 

±2.4xl0-2 

±1.3xl01 

±9.1x10-2 

Increasing  the  spatial  resolution  is  an  effective  way  to  reduce  the  uncertainty  in  the  PIV 
measurements.  For  example  improving  the  spatial  resolution  by  a  factor  of  two  (i.e.  an 
image  magnification  increase  of  50%)  will  decrease  the  instantaneous  axial  velocity 
uncertainty  by  approximately  40%.  A  higher  spatial  resolution  can  be  achieved  by 
mounting  the  same  camera  inside  the  test  section  of  the  wind  tunnel  or  by  using  a  lens 
with  a  longer  focal  length.  The  drawback  with  either  of  these  approaches  would  be  the 
reduced  field-of-view  and  therefore  in  order  to  achieve  the  same  viewing  area  more  field- 
of-views  would  be  needed  to  be  blended  together. 

The  uncertainty  in  the  mean  parameters  can  be  reduced  by  taking  more  samples,  N. 
Assuming  the  standard  deviation  in  the  velocity  measurements  remains  the  same,  taking 
3000  samples  instead  of  1000  would  reduce  the  uncertainty  in  the  mean  velocity,  vorticity 
and  Reynolds  stress  by  approximately  40%  (refer  to  Appendix  C). 

4.4.2  Convergence  of  Mean  Flow  Parameters 

The  convergence  of  the  mean  flow  was  assessed  qualitatively  by  plotting  a  running 
average  for  each  flow  parameter.  Plots  of  the  running  average  of  the  streamwise  velocity, 
u,  and  the  vertical  velocity,  v,  for  a  =  15°  and  field-of-view  3,  at  a  location  {x/L,  y/L}  = 
{-0.0280,  0.0025}  relative  to  the  tail  tip,  are  shown  in  Figure  25.  The  mean  velocities 
stabilise  and  approach  an  approximately  constant  value  after  -200  samples.  With 
increasing  sample  size,  the  mean  velocities  remain  statistically  stationary  in  time  with  no 
accumulative  bias  uncertainty  evident  in  the  mean.  The  normalised  ^-component  velocity 
fluctuation  ( u'/U^),  v-component  velocity  fluctuation  (v'/U^)  and  Reynolds  stress 

(  u'v' /U ^  )  running  averages  are  also  presented  in  Figure  25,  and  stabilise  after 
approximately  500  samples. 
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The  convergence  of  the  data  at  several  other  locations  in  the  flow  was  also  examined, 
including,  near  the  surface  of  the  model  and  in  the  wake  region.  Similar  behaviours  were 
observed,  with  convergence  of  all  parameters  occurring  within  approximately  500 
samples.  These  results  further  validate  that  the  number  of  image  fields  taken 
(approximately  1000)  for  each  test  case  was  adequate  in  establishing  convergence  in  the 
statistical  parameters. 


Number  of  fields 

Figure  25  -  Running  average  plots  of  (a)  streamwise  velocity  (u),  (b)  vertical  velocity  (v)  and  (c) 
normalised  u-component  velocity  fluctuation  (u1  /uj),  v-component  velocity  fluctuation 

(V/Uj  and  Reynolds  stress  (u'v'/UD  at  a=  15°  for field-of-view  3  at  fx/L,  y/L}  =  {- 
0.0280,  0.0025}. 

4.4.3  Assessment  of  Peak-Locking 

Following  the  procedure  in  Raffel  et  al.  [4],  for  the  image  capture  settings  used  in  this 
study  (i.e.  lens  f-stop  of  8  and  a  reproduction  ratio  of  6.67),  the  diffraction  limited  particle 
image  diameter  was  estimated  to  be  1.33  pixels.  Peak-locking  was  assessed  by  plotting 
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histograms  of  particle  displacement.  Figure  26  is  a  typical  histogram  plot  showing  the 
particle  displacement  counts  at  a  =  -2.5°.  The  displacement  counts  are  not  biased  towards 
integer  pixel  values  indicating  peak-locking  was  not  present.  The  histogram  plots  from  the 
other  test  cases  also  showed  an  acceptable  distribution. 


Figure  26  -  Histogram  ofPIV  displacement  data  at  a  =  -2.5° and  field-of -view  3,  indicating 
unbiased  pixel  displacement.  Histogram  bin-width  =  0.05  pixel. 


4.5  Lessons  Learnt  and  Key  Recommendations 

PIV  is  a  new  technique  being  employed  in  the  DSTO  experimental  facilities,  and  its 
implementation  in  a  large-scale  facility  such  as  the  LSWT  presents  a  number  of  technical 
challenges.  This  section  reports  the  lessons  learnt  during  these  tests,  with  a  view  towards 
providing  guidance  and  improvements  for  future  test  campaigns. 

4.5.1  Surface  Reflections 

The  high  laser  light  intensities  used  in  the  experiment  produce  strong  reflections  which,  in 
some  cases,  resulted  in  invalid  measurement  regions  close  to  the  surface  of  the  model.  The 
strong  reflections,  if  directed  back  towards  the  lens,  also  have  the  potential  to  destroy  the 
pixel  elements  in  the  CCD  camera.  Furthermore,  the  reflected  light  could  illuminate 
particles  outside  the  incident  light  sheet  resulting  in  spurious  measurements. 

The  submarine  model  was  not  specifically  designed  for  PIV.  The  aluminium  surface  was 
polished  and  then  anodised  resulting  in  a  highly  reflective  surface  finish.  Initial  tests  at 
low  laser  power  levels  revealed  that  significant  image  flaring  in  regions  near  the  model 
surface.  Furthermore,  due  to  the  curvature  of  the  tail,  reflections  towards  the  CCD  camera 
were  of  particular  concern.  In  order  to  avoid  direct  reflections  back  to  the  camera  the  laser 
sheet  was  shifted  in  the  starboard  direction  (positive  z-direction)  away  from  the  centre-line 
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of  the  model  by  approximately  1  mm.  This  reduced  the  image  flare-up  near  the  surface  of 
the  model.  For  future  PIV  test  campaigns,  the  model  surface  should  be  bead  blasted  before 
anodising  to  ensure  a  matt  finish  in  order  to  reduce  effects  of  surface  reflections. 

4.5.2  Optical  Quality  of  Test  Section  Window 

Initial  PIV  shakedown  tests  highlighted  the  poor  optical  quality  of  the  LSWT  working 
section  window.  This  proved  to  be  an  impediment  in  obtaining  consistent  high  quality  PIV 
data.  The  poor  optical  quality  of  the  window  produced  optical  distortions  and  are  most 
likely  due  to  an  inhomogeneous  material,  varying  thickness  of  the  glass,  and  surface 
curvature.  The  test  section  window  also  contained  a  number  of  scratches  and  marks. 
Taking  PIV  measurements  through  these  regions  of  poor  optical  quality  can  result  in  the 
olive  oil  seed  particles  appearing  out-of-focus.  This  reduces  the  resultant  cross-correlation 
peak-to-noise  ratio  and  produces  a  high  proportion  of  invalid  data  in  these  areas. 

Due  to  the  poor  optical  quality  of  the  LSWT  test  section  window  the  following  procedural 
changes  were  adopted: 

1.  The  camera  was  initially  positioned  at  a  desired  location  such  that  the  field-of-view 
corresponded  to  one  of  the  regions  of  interest  depicted  in  Figure  4. 

2.  An  initial  measurement  was  then  taken  and  the  quality  of  the  PIV  data  was 
assessed.  For  the  majority  of  test  cases  this  initial  measurement  produced 
unsatisfactory  PIV  data  with  large  areas  of  invalid  data  (e.g.  the  spurious  vector 
count  was  greater  than  5%). 

3.  The  camera  was  then  translated  on  the  traverse  to  a  new  location  and  the 
measurements  repeated.  The  camera  was  moved  until  the  window  was  of 
sufficient  optical  quality  to  produce  good  PIV  measurements  (i.e.  the  spurious 
vector  count  was  less  than  5%). 

For  some  cases  it  was  not  possible  to  obtain  satisfactory  PIV  measurements  over  the  entire 
field-of-view.  In  these  cases  the  regions  of  the  image  which  consistently  produced 
spurious  vectors  were  masked  out. 

In  future  PIV  test  campaigns  the  replacement  of  the  test  section  windows  using  a  material 
with  better  optical  grade  is  recommended.  Candidate  materials  include  float  glass  or 
aircraft  grade  Plexiglas;  the  exact  specification  and  quality  of  such  materials  and  their 
suitability  for  the  current  application  is  matter  for  further  investigation. 

4.5.3  Seeding  Losses 

The  DSTO  LSWT  is  open  to  atmospheric  pressure  and  therefore  seeding  losses  can  occur 
due  to  leakage  through  gaps  in  the  circuit.  Losses  also  occur  by  oil  depositing  on  exposed 
surfaces  such  as  the  model  surface,  support  structure,  tunnel  walls,  turning  vanes,  screen 
and  the  fan. 

Raffel  et  al.  [4]  showed  that  the  number  of  particles  in  each  IW  must  be  greater  than  8-12  to 
ensure  the  valid  detection  probability  exceeds  95%  and  therefore  these  losses  have  the 
potential  to  degrade  the  quality  of  the  PIV  data.  The  operation  of  the  seeding  system 
showed  direct  correlation  with  the  amount  of  spurious  vectors  in  the  PIV  results.  Figure  27 
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shows  two  plots  of  the  spurious  vector  count  over  the  number  of  fields  taken  at  (a)  a  =  0° 
and  (b)  a  =  5°.  At  a  =  0°,  the  seeding  system  was  deliberately  toggled  ON  and  OFF  and  as 
a  consequence  the  spurious  vector  count  shows  distinct  peaks  and  troughs.  At  a  =  5°,  the 
seeding  system  remained  ON  during  the  entire  run  and  as  a  consequence  the  spurious 
vector  count  remains  low.  For  the  test  conditions  (Uoo  =  30  m/  s)  and  the  seeding  system 
ON  during  the  entire  image  capture  sequence,  seeding  losses  were  accounted  for  and 
adequate  seeding  remained  inside  the  tunnel  without  over  saturation.  These  results  show 
that  in  order  to  obtain  consistently  good  quality  PIV  measurements,  for  a  range  of  test 
conditions,  it  is  necessary  for  the  seeding  concentration  to  be  monitored  such  that  the 
minimum  acceptable  level  can  be  maintained. 

After  several  days  of  testing  the  test  section  walls  and  the  model  surface  were  covered 
with  a  film  of  oil  residue.  Testing  at  wind  tunnel  velocities  greater  than  30  m/  s  resulted  in 
increased  rate  of  oil  deposition  on  the  tunnel  surfaces  resulting  in  significant  seeding 
losses.  Seeding  improvements  may  be  achieved  via  local  seeding  injection  or  an  upstream 
rake.  Both  of  these  methods  are  intrusive  and  their  effects  on  the  flow-field  would  need  to 
be  assessed.  In  order  to  avoid  oil  build  up  on  exposed  surfaces,  other  seeding  alternatives 
should  be  investigated.  One  alternative  to  olive  oil  seeding  is  Di-Ethyl-Hexyl-Sebacate 
(DEHS)  which  is  a  non-soluble  liquid  suitable  for  producing  aerosols  and  has  been  used 
for  PIV  seeding  in  other  international  large-scale  facilities  [4]  and  subject  of  further 
investigation. 
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Figure  27  -  Plot  of  spurious  vector  count  over  the  number  of  fields  taken  at  (a)  a=  0°  and  (b) 
a  =  5°,  for  field-of-view  1  at  afreestream  airspeed  of  30  m/s. 
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5.  Conclusions 


DSTO  researchers  have  utilised  the  complimentary  techniques  of  smoke  flow  visualisation 
and  PIV  to  investigate  the  flow-field  over  a  generic  submarine  model  tested  in  the  LSWT. 
Smoke  flow  visualisation  revealed  the  large-scale  structure  of  the  flow-field  around  the 
model.  The  flow-field  was  characterised  by  cross-flow  separation  that  forms  near  the  tail, 
and  migrated  forward  with  increasing  incidence. 

The  PIV  technique  provided  mean  flow-fields,  vorticity  maps  and  contour  plots  of 
Reynolds  stress  at  targeted  areas  of  interest.  With  increasing  incidence,  the  streamlines 
over  the  upper  surface  were  vectored  down  along  the  surface.  The  PIV  results  showed 
flow  features  consistent  with  the  flow  visualisation  results.  The  PIV  vorticity  contours 
showed  that  with  increasing  incidence,  the  boundary  layer  vorticity  on  the  upper  surface 
extended  further  downstream. 

Both  the  flow  visualisation  and  PIV  measurements  showed  that  the  model  mounting 
arrangement  produced  flow  disturbance  on  the  underside  of  the  model.  Despite  this 
interference  the  large-scale  flow  features  were  consistent  with  those  reported  for  the  flow 
over  a  prolate  spheroid. 

The  PIV  results  were  compared  with  CFD  results  and  a  reasonable  agreement  in  the 
streamwise  velocity  and  vorticity  were  obtained.  Given  that  both  techniques  produced 
similar  results,  the  comparison  with  CFD  was  both  informative  and  served  to  validate  the 
experimental  method. 

The  PIV  testing  provided  insight  into  the  complexities  involved  in  obtaining  PIV 
measurements  in  a  large-scale  wind  tunnel  facility.  Practical  limitations  of  the  current  PIV 
technique  and  the  facility  dictated  the  test  conditions  for  which  good  measurements  could 
be  acquired.  The  optical  quality  of  the  test  section  window  was  highlighted  as  a  deficiency 
in  the  experimental  setup,  as  was  the  importance  of  consistent  seeding,  in  obtaining  high 
quality  PIV  results. 

A  range  of  improvements  have  been  suggested  which  are  currently  being  implemented. 
With  these  refinements  applied  to  the  PIV  test  technique,  researches  are  confident  that  an 
enhanced  PIV  capability  will  be  provided  for  future  test  campaigns. 
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Appendix  A:  Wind  Speed  Measurement 

The  DSTO  LSWT  flow  reference  instrumentation  was  used  for  measurement  of  the  tunnel 
static,  the  differential  pressure  across  the  contraction  and  the  tunnel  temperature.  These 
measurements  are  used  to  calculate  the  test  section  wind  speed  and  Reynolds  number. 

The  differential  pressure  is  obtained  by  measurement  of  the  pressure  between  two 
piezometer  rings  located  at  the  entrance  and  exit  of  the  tunnel  contraction.  The  pressure  is 
sampled  by  a  3  PSID  range  Digiquartz  differential  pressure  transducer.  A  calibration  factor 
(1.079)  is  applied  as  a  correction  to  the  pressure  measurement  obtained  from  the  two 
piezometer  rings  to  represent  accurately  the  correct  dynamic  pressure  and  consequently, 
velocity,  at  the  centre  of  the  test  section  midway  between  the  centres  of  the  turntables  in 
the  floor  and  ceiling  [23].  The  static  pressure  in  the  test  section  is  measured  by  an  absolute 
23  PSIA  Digiquartz  pressure  transducer.  The  tunnel  temperature  is  measured  by  an  Intech 
LPN-R  RTD  temperature  transmitter  in  combination  with  a  3  wire  Platinum  PtlOO  RTD 
temperature  transducer. 


A.l.  Wind  Speed  Measurement  Uncertainty 

The  measured  parameters  required  to  estimate  the  wind  speed  in  the  test  section  of  the 
LSWT  are 


Pd  =  Differential  pressure  across  the  piezoelectric  rings  (Pa), 
Ps  =  Static  pressure  of  the  test  section  (Pa), 

T  =  Test  section  temperature  (K). 

Other  parameters  required  in  the  calculation  of  the  wind  speed  are 


P  =  Density  of  air  in  the  test  section  (kg.nv3), 

k  =  Wind  speed  calibration  factor  of  the  tunnel  =  1.079  [23], 

R  =  Ideal  gas  constant  for  air  =  287  (J.kg^K-1). 


The  wind  speed  or  freestream  airspeed,  U  f  is  calculated  according  to  the  equation. 


U„ 


(A.1) 


Where  the  air  density,  p ,  is  evaluated  according  to  the  ideal  gas  equation. 


P  = 


RT' 


(A.2) 
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The  uncertainty  in  U ^  is  then  estimated  by  combining  the  differential  pressure 
uncertainty,  Epdf  the  wind  speed  calibration  factor  uncertainty,  E*>  and  the  uncertainty  in 
the  density,  Ep , 


or,  when  simplified. 


dU  oo 
KdPd 


JPd 


+ 


dlL 

dk 


2  ( 
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dU, 


Y 


dp  Ep 


EU  = 


fkEl' 1  f  PdE\  ^  A 


2p^ 
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2p& 


+ 


kPdE2p  ^ 


2p3 


(A.3) 


(A.4) 


where  the  uncertainty  in  the  density  is  estimated  by  combining  the  statics  pressure 
uncertainty,  Eps  and  the  tunnel  temperature  uncertainty,  Ep , 


f  ~  'N 
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l  RT 2  T) 

(A.5) 


For  typical  tunnel  conditions,  the  uncertainties  based  on  random  variable  analysis  for  a 
95%  confidence  interval,  that  is  a  uncertainty  coverage  factor  of  2,  are  as  follows. 


101325  Pa, 
300  K, 


=  1.079, 


Eps 

Et 


±  6.10  Pa  [24], 

±  0.5  K  (based  on  a  conservative  accuracy 
estimate  for  a  resistance  temperature  detector). 


±0.017,  U  <20 m/s, 


[23] 


±0.003,  Ux>40m/s, 

For  30  m/s,  an  average  of  Ek  =  ±  0.01  is  assumed. 


Epd 


2.18  Pa  [25], 


For  the  typical  tunnel  conditions,  the  wind  speed  uncertainties  for  a  range  of  tunnel 
velocities  are  shown  in  Table  Al. 


Table  Al  -  LSWT  wind  speed  uncertainties 


Wind  Speed 
llo,  (m/s) 

Differential  Pressure 

Pd  (Pa) 

Uncertainty 

Eu oo  (m/s) 

5.0 

13.6 

±0.40 

20.0 

218.1 

±0.19 

30.0 

490.8 

±0.16 

40.0 

872.5 

±0.08 
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Appendix  B:  PIV  Measurement  Settings 


Table  B1  -  Summary  of  PIV  measurement  settings 


Laser  Maximum  Pulse  Energy 

200  mj 

Laser  Wavelength 

532  nm 

Laser  Sheet  Thickness  (at  the  model  surface) 

1.5  mm 

Camera  Resolution 

4008  x  2672  pixels 

Camera  Pixel  Size 

9  jLim  x  9  jam 

Camera  Frame  Rate 

2.07  Hz 

Time  Delay  (At) 

6  - 15  ps 

Lens  Focal  Length 

200  mm 

Spatial  Resolution 

60.03  pm/ pixel 

First  Interrogation  Window  Size 

64  x  64  pixels 

Second  Interrogation  Window  Size 

64  x  64  pixels 

Overlap  between  Interrogation  Windows 

50% 

Vector  Spatial  Resolution 

1.92  mm 

Seed  Particles 

Olive  Oil 

Nominal  Particle  Size 

1  pm 

Number  of  Fields  (N) 

-1000 
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Appendix  C:  PIV  Uncertainty  Analysis 


C.l.  PIV  Measurement  Uncertainty 


The  PIV  measurement  uncertainties  in  the  instantaneous,  mean  and  fluctuating  velocity 
components  were  calculated  using  the  approach  described  in  [21].  The  total  PIV  velocity 
measurement  uncertainty,  Eu ,  is  obtained  by  combining  the  PIV  particle  displacement 
uncertainty  (Edisp),  scaling  uncertainty  ( ESCaie )  and  timing  uncertainty  (Etime), 


E2  = 


?2 

J  disp 


?2 

J scale 


+  Etime  ~ 


AX 


scale  J 


vA  tj 


(Cl) 


The  theoretical  minimum  resolvable  pixel  displacement,  6ax,  for  the  PIV  algorithm  was 
determined  to  be  0.05  pixel  for  the  recursive  algorithm.  This  was  obtained  by  processing 
digitally  simulated  images  with  known  particle  displacement.  The  values  quoted  in  the 
literature  are  of  the  same  order,  for  example,  a  minimum  resolvable  displacement  of 
0.05  pixel  was  also  attained  in  [22].  For  the  image  capture  settings  used  in  this  study,  a 
freestream  velocity  of  30  m/s  and  a  time  delay.  At  =  10  ps  corresponds  to  a  displacement 
of  5  pixels  (AX),  and  therefore  the  PIV  displacement  uncertainty  expressed  as  a  percentage 
of  the  freestream  velocity  is  ±0.5%  (1.0%). 

The  PIV  camera  images  are  calibrated  by  inserting  a  ruler  in-line  with  the  laser  sheet  and 
measuring  the  number  of  pixels  over  a  calibration  length.  For  the  image  calibration,  the  in¬ 
plane  scaling  uncertainty  was  estimated  to  be  2  pixels  based  on  the  ability  to  estimate  the 
centre  of  the  reference  graduations  on  a  ruler  from  an  image  acquired  using  the  PIV 
camera.  Two  graduations  need  to  be  chosen  over  a  calibration  length  thus  for  each 
graduations  there  is  2  pixel  uncertainty  in  both  the  horizontal  and  vertical  directions. 
Therefore  the  combined  total  in-plane  scaling  uncertainty  is 

Sscale  =  V(22  +22)  +  (22  +22)  =  4  pixels. 


For  the  out-of-plane  scaling  uncertainty  an  Empire™  spirit  level  was  used  to  align  the 
calibration  ruler  in  both  pitch  and  roll.  The  accuracy  of  the  spirit  level  is  quoted  as 
0.5  mm/ m  thus  for  the  6" (152.4  mm)  spirit  level  this  corresponds  to  0.0762  mm  (or  0.029°). 
The  projected  out-of-plane  uncertainty  is  then  equal  to 

152.4  -  152.4  cos  (0.029°)  =  1.92  xl0'5mm, 

and  given  the  magnitude  of  this  uncertainty  it  was  deemed  insignificant  and  was  ignored. 

The  calibration  length  used  was  120  mm.  For  an  image  spatial  resolution  of  60.0  pm/pixel 
this  corresponds  to  2000  pixels  (Lscaie).  Thus  the  total  scaling  uncertainty  is, 

Escaie  =  6scaie  /LSCaie  =  4/2000  =  0.002  =  ±0.1%  (total  uncertainty  =  0.2%). 
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This  analysis  assumes  the  ruler  is  calibrated  and  does  not  take  into  account  the  error  of  the 
actual  ruler  measurement  length. 

The  timing  uncertainty  is  a  combination  of  the  minimum  clock  time  of  the  synchroniser 
(1  ns)  and  the  laser  jitter  (1  ns).  The  combined  timing  uncertainty  for  two  laser  pulses  can 
be  calculated. 

For  a  time  delay.  At  =  10  ps  =  10000  ns,  the  corresponding  timing  uncertainty, 

Ew  =  6m  /  At  =  2/10000  =  ±  0.0002  =  ±  0.02%  (total  uncertainty  =  0.04%). 


The  total  PIV  velocity  uncertainty  for  the  u  velocity  component  is  then  given  by  the  root- 
sum-square  of  each  the  uncertainties, 

Eu  =JeIp  +£ Lle  +^L  =  a/0.52  +  0. 1 02  +0.022  ~±0.51%(\.02%) . 

Therefore  the  normalised  instantaneous  PIV  streamwise  velocity  (u/Uk)  the  corresponding 
uncertainty  is  ±0.0051. 

The  combined  PIV  streamwise  velocity  and  LSWT  freestream  airspeed  uncertainty  is 
found  by  a  summation  in  quadrature  of  the  corresponding  uncertainties. 


( Eu ± 

2 

,2  f  0.16  V 

£— =  1 

\e;  + 

=  ^ 

,  0.00512  + 

i 

l  Vo j 

V  l  30  J 

=  ±0.74%  (T. 48%/ 


Therefore  for  the  normalised  instantaneous  PIV  streamwise  velocity  ( u/Uoo )  the  combined 
uncertainty  is  ±0.0074. 

The  v  velocity  component  is  assumed  to  have  identical  uncertainty  as  the  u  velocity 
component,  i.e.  Ev  -  ±0.0051  (0.0102)  andi?^-  =  ±0.0074  (0.0148). 


C.2.  Mean  Velocity  Uncertainty 

Velocity  measurements  obtained  from  PIV  are  considered  to  be  statistically  independent 
and  conform  to  a  normal  distribution  [22].  From  the  theory  of  random  variable  analysis, 
for  a  95%  confidence  interval,  the  fractional  uncertainty  in  the  streamwise  (u)  mean 
velocity  (or  the  statistical  uncertainty)  is  given  by. 


E-  = 


2<g»  > 

VoW 


9  u 

RMS 

U„y[N’ 


(C.2) 


44 


UNCLASSIFIED 


UNCLASSIFIED 


DSTO-TR-2944 


where  ou  is  the  standard  deviation,  urms/Ux,  is  the  turbulence  intensity  and  N  is  the 
number  of  samples.  The  existence  of  a  PIV  velocity  uncertainty,  Eu,  in  each  individual 
measurement  is  added  to  the  statistical  uncertainty,  Eu,  to  get  the  total  uncertainty  in  the 
mean  velocity.  The  maximum  turbulence  intensity  (uRMs/Ux)max  was  determined  from  the 
mean  flow-fields  and  was  found  to  be  0.34,  therefore 

_  2x0,34 

“  Viooo 

=  ±2.2%  (4.4%) 

The  total  uncertainty  in  the  streamwise  (u)  mean  velocity  is  then  given  by  the  root-sum- 
square, 

E  -  =  J(E-  7  +  El  =  J(  2.2/  +0.512 

Total  u  V  u  '  u  V '  ' 

-  ±2.3%  (4.5%). 


Therefore  for  the  normalised  mean  PIV  derived  velocity  ( U/U o0)  the  corresponding 
uncertainty  is  ±0.023. 

The  combined  normalised  mean  PIV  velocity  and  freestream  airspeed  uncertainty  is  found 
by  a  summation  in  quadrature  of  the  corresponding  uncertainties. 


E 


u 


E2  ,-  + 

( E  \ 

Total  u 

u 

V  00  J 

±2.4%.  (4.8%) 


^0.023 2 


f  0,16' 
l~30~. 


Thus  for  the  normalised  mean  PIV  velocity  (u/Uf  ) the  corresponding  combined  fractional 
uncertainty  is  ±0.024. 

The  mean  v  velocity  component  is  assumed  to  have  identical  uncertainty  as  the  u  mean 
velocity  component. 


C.3.  Vorticity  Uncertainty 

Two  uncertainties  make  up  the  out-of-plane  vorticity  uncertainty,  the  precision  (or 
random)  uncertainty  and  the  bias  uncertainty  [7],  The  precision  uncertainty  represents  the 
uncertainty  due  to  velocity  field  calculation  (i.e.  the  PIV  uncertainty)  and  the  bias 
uncertainty  represents  the  consistent  under-estimation  or  over-estimation  of  the  vorticity 
estimation  scheme.  Both  depend  on  the  spatial  separation  of  the  data.  A,  and  on  the  size  of 
a  relevant  length  scale  to  be  resolved,  Lc,  although  in  an  opposing  fashion. 

The  normalised  precision  uncertainty  for  the  second  order  finite  difference  vorticity 
calculation  scheme  used  in  this  study  may  be  calculated  from  [7]  as, 

EVes=^Eu,  (C.3) 
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where  Eu  is  the  total  PIV  velocity  uncertainty  in  the  streamwise  velocity  and  A  is  a  non- 
dimensional  precision  uncertainty  transmission  ratio  and  is  given  as 

A  =  k^,  (C.4) 

A 

where  k  =  1  for  the  second  order  finite  difference  scheme.  In  this  analysis  it  is  assumed 
that  the  v  velocity  component  uncertainty  is  the  same  as  the  u  velocity  uncertainty. 

The  normalised  peak  vorticity  bias  uncertainty  for  the  second  order  finite  difference 
scheme  can  be  calculated  using  [7], 


(C.5) 


where  oomax  is  the  peak  out-of-plane  vorticity  calculated  from  the  PIV  velocity 
measurements  and  oobias  is  the  actual  peak  out-of-plane  vorticity  value,  y  =  0.2416  and 
(/)  =  2  for  the  second  order  finite  difference  scheme. 


Thus  the  total  uncertainty  in  the  vorticity  is  calculated  by. 


(C.6) 


Where  the  total  uncertainty  in  the  vorticity  is  expressed  as  fractional  uncertainty  of  the 
actual  out-of-plane  vorticity  and  the  peak  out-of-plane  vorticity  calculated  form  the  PIV 
velocity  measurement,  i.e.. 


Substituting  Eq.  (C.3)  and  Eq.  (C.5)  gives, 

(C.7) 

For  a  64  x  64  IW  size,  50%  overlap  factor  and  an  image  spatial  resolution  of  60.03  pm,  the 
spatial  separation  of  the  data  is  A  =  32  x  0.06003  =  1.92  mm.  The  total  PIV  uncertainty  in 
the  streamwise  velocity  was  calculated  previously,  Eu«  ±0.51%. 

As  a  representative  case  to  explore  the  effect  of  spatial  resolution  and  the  vorticity  estimate 
scheme  on  the  bias  and  random  error,  Etabari  and  Vlachos  [7]  use  an  Oseen  vortex  and 
defined  Lc  as  the  as  length  scale  of  the  vortex  core.  The  flow  around  submarine  can  be 
considered  a  complex  flow  field  with  a  range  of  length  scales  present  in  the  flow,  thus  a 
characteristics  length  scale  is  difficult  to  define.  A  conservative  assumption  was  made  and 
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the  characteristic  length  scale  was  based  on  the  vorticity  thickness  at  the  tail  of  the  model, 
which  is  defined  as 


U  MAX  U  min 


r  dll' 

(  du' 

’ 

(C.8) 

[dy) 

MAX 

Kdy) 

MAX 

where  Umin  =  0  at  the  surface  of  the  model.  The  vorticity  thickness  at  the  tail  of  the  model 
was  determined  to  be  approximately  21  mm. 

Thus  the  uncertainty  in  the  instantaneous  vorticity  is 


E 


CO 


21 

L92 


\2  (  /, 


0.0051 


j 


0.2416 


V 


1.92 

~2\ 


±0.056 
±5. 6%  (11%) 


The  fractional  uncertainty  in  the  instantaneous  vorticity  is  therefore  +0.056  of  the 
maximum  PIV  measured  instantaneous  vorticity. 

The  mean  vorticity  components  were  extracted  from  the  mean  velocity  field  and  the 
uncertainty  calculation  method  is  thus  identical  to  that  explained  above.  The  bias 
uncertainty  remains  the  same  while  the  precision  uncertainty  increases  in  accordance  with 
the  greater  uncertainty  associated  with  the  mean  velocity  field,  i.e. 

E^„=A-iE:+El)  <C9> 


The  value  A  is  the  same  as  (C.4)  and  the  mean  velocity  uncertainty  in  the  streamwise 
velocity  was  previously  calculated,  E-  =  +2.3%,  thus  the  total  uncertainty  in  the  mean  out- 
of-plane  vorticity  is  calculated  using. 


E- 
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XfiTX) )  + 


f 


0.2416 


'a' 


v 


\Lc  j 


Substituting  the  values. 


(C.10) 
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1.92 
±0.25 

±25%  (50%). 
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The  fractional  uncertainty  in  the  mean  out-of-plane  vorticity  is  calculated  to  be  ±0.25  of  the 
maximum  measured  PIV  vorticity. 

The  combined  normalised  mean  PIV  vorticity,  model  length  and  freestream  airspeed 
uncertainty  is  found  by  a  summation  in  quadrature  of  the  corresponding  uncertainties. 


f  17  A 


E =JE -  + 


c oL 

El 


+ 


V  ^  J 


=  J0.252  + 


0.5 

1350 


+ 


0.16 
30  j 


=  ±0.25 
=  ±25%  (50%). 


where  the  uncertainty  in  the  model  length,  El  =  ±0.5  mm. 

The  maximum  normalised  mean  PIV  vorticity,  (aZMAXL/UJ  *  53  and  thus  the 
corresponding  mean  out-of-plane  vorticity  uncertainty  is  E—  ~  (0.25  x  53)  «  ±13. 

ull 


C.4.  Reynolds  Stress  Uncertainty 

The  normal  stresses,  u'2  and  v'2  are  related  to  the  standard  deviation  of  each  velocity 
components  as  follows. 


u 


,2 


V 


|2 


(C.U) 

(C.12) 


and  from  the  theory  of  propagation  of  uncertainties  [26],  their  associated  fractional 
uncertainties  are 


^=2  E„u, 


E~  =  2E„  ■ 

i2  (7  V 


(C.13) 

(C.14) 


For  a  95%  confidence  interval  the  corresponding  fractional  uncertainty  in  the  standard 
deviation  of  the  velocity  may  be  written  [27], 


4W2 1)’ 

2 

Mn-i) 


(C.15) 

(C.16) 
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The  values  of  Eou  and  Eov  are  thus  independent  of  the  sample's  standard  deviation  and 
merely  a  function  of  sample  size  [27].  The  Reynolds  shear  stress  is  considered  to  have 
similar  uncertainties  associated  with  them  and  the  uncertainties  associated  with  the 
velocity  components  can  be  considered  equal,  i.e.  E- 2  =  E_2  =  E— ,  therefore, 

u'  v'  u  v 


E—  = 

U  V 


4 

J2(N-1)' 


(C.16) 


Thus  for  a  sample  size  of  N  =  1000,  the  uncertainty  in  the  ensemble-averaged  Reynolds 
stress  (and  normal  stresses)  is  ±9%  (total  uncertainty  =  18%)  as  a  percentage  of  the 
measured  value.  Thus  for  the  fractional  uncertainty  of  the  PIV  derived  normalised 


Reynolds  stress. 


u  v 


KUlj 


,  the  corresponding  uncertainty  is  ±0.09. 


The  combined  normalised  PIV  Reynolds  Stress  and  freestream  airspeed  uncertainty  is 
found  by. 
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u  v 
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(e—)-  + 

V  u  v  / 


2  E, 


V  U«  7 


0.092  + 


^  2x016  V 

v  30  , 


±0.091 

±9.1%  (18.2%). 


Thus  for  the  normalised  Reynolds  stress. 


mV 
U 2 


,  the  corresponding  combined  fractional 


uncertainty  is  ±0.091. 
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